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THRUST COEFFICIENTS OF LOW-THRUST NOZZLES 


by Ernie W. Spisz, Paul F. Brinich, and John R. Jack 
Lewis Research Center 
SUMMARY 

The performance losses associated with small nozzles operating at low propellant 

flow rates were investigated experimentally with a resistance-heated hydrogen thrustor. 

Data were obtained for seven different nozzles covering the range of nozzle area ratio 

-4 -4 -4 

from 25 to 150, hydrogen propellant flow rates of 0. 5x10 , 1. 0x10 , and 2. 5x10 

pound per second, and propellant temperatures from 530° to 4000° R. 

The nozzle losses, as determined by the difference between the one-dimensional 
isentropic thrust coefficient and the measured thrust coefficient, were correlated in 
terms of the throat Reynolds number and nozzle area ratio. At low throat Reynolds num- 
ber, the nozzle losses were sufficiently large to result in measured thrust coefficients 
less than the values calculated for isentropic flow through a choked orifice. Results 
show that large area ratios are not required for achieving maximum thrust for nozzles 
with low throat Reynolds number (<10 000). 


INTRODUCTION 

It has been generally recognized (refs. 1 and 2) that viscous boundary-layer and heat- 
transfer effects on the divergent surfaces of supersonic nozzles can degrade the per- 
formance (thrust, specific impulse, and efficiency) of low-thrust thermal-propulsion de- 
vices. The performance losses associated with these nozzle -expansion effects have 
not been of immediate concern because they are negligible for the relatively high thrust 
levels (>500 millipounds) of thrustors that are being developed. The experimental per- 
formance values of these thrustors are generally within a few percent of the calculated 
values for a one-dimensional isentropic nozzle-expansion process. (Typical one- 
dimensional isentropic nozzle-expansion calculations for hydrogen are presented in 
refs. 3 and 4.) However, for lower thrust levels (<100 millipounds, which may be de- 
sirable for applications such as attitude control of space vehicles) the boundary-layer and 



heat-transfer effects that occur in the divergent portion of the nozzle become more pro- 
nounced. For such small, low-thrust nozzles, the boundary layer can account for a sub- 
stantial portion of the available nozzle flow area. The boundary layer introduces viscous 
losses and alters the free-stream conditions so that the conventional assumption of a one- 
dimensional isentropic expansion process becomes questionable. The nozzle flow is fur- 
ther modified as a result of the thermal boundary layer and the corresponding heat- 
transfer effects that are intimately coupled to both the free-stream and boundary-layer 
flows. For low propellant flow rates, heat-transfer losses to the nozzle walls can be a 
significant portion of the total thermal energy available for conversion to kinetic energy. 

Little experimental data are available for estimating the boundary-layer and heat- 
transfer losses at low-thrust levels. Some overall performance data are presented in 
references 5 and 6, but these data primarily indicate only the performance losses that 
can occur at low-thrust levels. The program herein was undertaken to study further the 
nozzle losses from performance data obtained with the radiation-cooled resistojet of ref- 
erence 7. Previous experience with this radiation-cooled resistojet at low-thrust levels 
had indicated that, as the propellant flow rate was decreased corresponding to a reduction 
in thrust level, the thrust coefficient also decreased. This decrease in thrust coefficient 
was interpreted in terms of a Reynolds number effect. This effect is investigated herein 
and covers nozzle throat Reynolds numbers ranging approximately from 500 to 18 000. 
This range was obtained by varying propellant flow rate, nozzle throat diameter, and vis- 
cosity (by propellant temperature changes). The effects of nozzle area ratio were studiecj 
for seven different nozzles with area ratios ranging from 25 to 150. All nozzle divergent 
sections were conical and had a divergence angle of 20°. Hydrogen was used exclusively 
as the propellant, and data were obtained for propellant temperatures ranging from 530° 
to 4000° R. 


PERFORMANCE CONSIDERATIONS 


For a conical nozzle such as that shown in figure 1, the thrust that would be de- 
veloped in a vacuum environment is given by 
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(All symbols are defined in appendix A. ) 
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Nozzle exit 

Figure 1. - Schematic drawing of nozzle. 



For a one -dimensional isentropic expansion process, equation (1) can be integrated 
and the thrust given by (ref. 8) 


F, = — u„A + p A. 
1 „ e *e e 

g 


( 2 ) 


where A is the nozzle angle divergence coefficient, A = i(l + cos a). 

If the thrust is nondimensionalized to obtain the thrust coefficient Cp, = F/P q A th’ 
the following equation for the thrust coefficient, in terms of the nozzle area ratio, pro- 
pellant specific heat ratio, and nozzle divergence angle, results: 
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To calculate the thrust for a nonisentropic and non-one-dimensional nozzle-expansion 
process, equation (1) requires that the distribution of density, velocity, and pressure at 
the nozzle exit be known. If this information is unavailable, it is necessary to resort to 
approximate calculation techniques. A technique often used is to define the thrust as the 
difference between the calculated isentropic thrust and an estimated thrust-loss term: 


F = Fj - F £ (4) 

For small nozzles, the thrust-loss term is composed primarily of viscous drag on the di- 
vergent surfaces of the nozzle and heat-transfer losses incurred during the expansion 
process: 


F l = 


F + 
v 


ht 


(5) 
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The thrust loss due to heat-transfer losses to the nozzle walls can be estimated 
for a specific nozzle geometry, cooling design, and operating condition. This loss, how- 
ever, is minor for radiation-cooled thrustors except at very low propellant flow rates and 
therefore will not be evaluated here. The thrust loss due to viscous effects F y is given 
by the integral of the local shear stress r over the divergent surfaces of the nozzle: 




2irrr cos a ds 


( 6 ) 


The magnitude of this loss can be estimated by assuming that the local shear stress 
is given by the skin friction relation for laminar flow over a flat plate (ref. 9). Although 
the laminar flat plate relation does not apply directly to axisymmetric nozzle flow with 
its inherent axial variations in the flow conditions and possible turbulence, the resulting 
estimates provide a simple first approximation for the comparison with experiment. 

With this assumption, the loss in thrust coefficient due to viscous drag on the sur- 
face of a conical nozzle was derived in appendix B and is given by 
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Figure 2. - Viscous losses as given by laminar-flow flat-plate theory. Stagnation temperature, 
530° R; wall- to stagnation-temperature ratio, 1.0. 


In order to evaluate C„ it is necessary to specify the wall temperature and the 

r v 

nozzle stagnation conditions. The local flow conditions throughout the nozzle can be de- 
termined approximately in terms of the stagnation conditions from one-dimensional 
isentropic-flow relations. The local stream conditions, M and T/T , in the factor fg 
may be evaluated at either the geometric nozzle area ratio (uncoupled case) or an ef- 
fective nozzle area ratio (coupled case). The coupled case attempts to include the in- 
fluence of the boundary layer on the local stream properties in terms of a reduced nozzle 
area ratio. For the coupled case, the effective nozzle area ratio is defined in terms of 
the displacement thickness 5* by 



where 5 */r is determined as outlined in appendix B. 

Figure 2 presents the results for the numerical integration of equation (7) for hydro- 
gen (y = 1. 4) at cold-flow conditions (T Q = 530° R, T w /T q « 1. 0) for both the uncoupled 
and coupled cases. As the value of the parameter (Re^ tan a) decreases, the loss in 
thrust coefficient due to viscous effects increases. 

Figure 3 presents the loss in thrust coefficient for a nozzle with a divergence angle 
of 20° for the coupled solution. The increasing losses resulting from increasing area 
ratio and decreasing throat Reynolds number are evident. 

Figure 4 presents the net thrust coefficient, including the viscous losses, as given in 
figure 3, for a nozzle divergence angle of 20°. This divergence angle was used because 
preliminary calculations showed that, for a given nozzle area ratio and throat Reynolds 
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Figure 3. - Loss in thrust coefficient due to viscous effects as given by 
laminar-flow flat-plate theory. Stagnation temperature, 530° R; wall- 
to stagnation-temperature ratio, 1.0; nozzle divergence angle, 20°. 



Figured. - Thrust coefficient including viscous losses as given by laminar-flow flat- 
plate theory. Stagnation temperature, 530° R; wall- to stagnation-temperature 
ratio, 1.0; nozzle divergence angle, 20°. 


number, the thrust coefficient is not greatly affected by changes in the divergence angle 
for 15 < of < 25°. The optimization of the thrust coefficient with area ratio, as shown in 
figure 4, occurs because the gains due to the increased area ratio are less than the addi- 
tional viscous losses imposed by the increased surface area. The significant loss in 
thrust coefficient and the low optimum area ratios at low throat Reynolds numbers are 
important factors that must be considered for the application of low-thrust devices. 
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EXPERIMENTAL APPARATUS 

The objective of this experimental program was to obtain performance data on a 
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Figure 5. - Wire coil resistance-heated thrustor. 


resistojet with hydrogen as the propellant at low-thrust levels for various nozzle area 
ratios and throat Reynolds numbers. The resistojet is reasonably well suited for provid- 
ing the necessary data in that its operating characteristics are inherently stable and re- 
liable over a wide range of operating conditions. A specified operating point can be set 
and maintained until equilibrium conditions are achieved. 


Thrustor Description 

The performance of resistance-heated hydrogen thrustors has been under investiga- 
tion for a number of years. Results of some of these studies are presented in references 
7, 10, and 11. Figure 5 illustrates the thrustor design used in reference 7 and in the 
present investigation. Basically, the design consists of a tungsten-wire coil heater in- 
serted into a nozzle body. The outer surface of the nozzle body is surrounded by a thick 
layer of insulation to minimize heat losses. The propellant is heated while passing di- 
rectly over the heater coil and is expelled through the nozzle. Table I indicates the seven 
different nozzle geometries used during the experimental program to provide the desired 
variations in nozzle area ratio and throat diameter. Two basic thrustor designs were 
used that differed only in size and method of supporting the heater coil. These two 
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TABLE I. - NOZZLE GEOMETRIES AND COLD- FLOW DATA 


Unit 

designa- 

tion 

(a) 

Nozzle 

area 

ratio, 

A e' /A th 

Nozzle 

throat 

flow 

area, 

A th’ 
sq in. 

Isentropic 

thrust 

coefficient, 

c *i 

Propellant 
flow rate, 
w, 

lb/sec 

Measured 

total 

pressure, 

P m> 

nr 

psia 

Measured 

thrust, 

F m> 

lb 

Measured 

thrust 

coefficient, 

C F 

m 

Difference 
between 
calculated 
isentropic 
and mea- 
sured 
thrust 
coefficient, 

C F 

F H 

Throat 

Reynolds 

number, 

Re th 

3W25 

25 

3. 63x10" 3 

1. 652 

2. 5x10" 4 

11. 5 

60. 9xl0" 3 

1.455 

0. 197 

11 000 





1.0 

4. 79 

22.9 

1. 318 

. 334 

4 400 





.5 

2. 44 

10. 57 

1. 192 

.460 

2 200 

3BN50 

50 

3. 84x10" 3 

1. 677 

2. 5X10" 4 

00 

tH 

61. 9xl0' 3 

1.442 

0. 235 

10 700 





1.0 

4. 60 

23.9 

1.353 

. 324 

4 280 





. 5 

2. 30 

10.78 

1. 222 

.455 

2 140 

3BN75 

75 

3. 63x10" 3 

1. 691 

2. 5X10" 4 

11.73 

61. 2x10" 3 

1.438 

0. 253 

11 000 





1.0 

4. 78 

23. 1 

1.331 

. 360 

4 400 





. 5 

2.48 

10. 58 

1. 175 

. 516 

2 200 

3BN100 

100 

4. 07X10' 3 

1. 696 

2. 5xl0“ 4 

10.3 

59. 8X10" 3 

1.425 

0. 271 

10 420 





1.0 

4. 39 

22.05 

1. 234 

.462 

4 170 





. 5 

2. 28 

10.02 

1.081 

. 615 

2 085 

4BN100 

100 

4. 65xl0" 3 

1. 696 

2. 5xl0~ 4 

8. 7 

61.7X10" 3 

1. 526 

0. 170 

9 730 





1.0 

3. 6 

22.7 

1. 355 

. 341 

3 890 





. 5 

1.9 

10.46 

1. 182 

. 514 

1 945 

4BN40 

40 

11.90X10" 3 

1. 670 

2. 5X10“ 4 

3. 55 

59. 8x10" 3 

1.415 

0. 255 

6 080 





1.0 

1. 50 

21. 85 

1. 224 

.446 

2 430 





. 5 

. 80 

10.02 

1.052 

.618 

1 215 

3BN150 

150 

1. 32x10" 3 

1. 705 

2. 5xl0” 4 

31.05 

59. 8x10" 3 

1.456 

0. 249 

18 330 





1.0 

12. 62 

22. 5 

1.350 

. 355 

7 320 





. 5 

6.48 

10. 14 

1. 188 

. 517 

3 665 


a First digit: 3 has 1/2- in. -diam heater coil and boron nitride heater supports, while 4 has 1-in. -diam heater 
coil and tungsten heater supports; letters W and BN indicate tungsten or boron nitride nozzle body; last 
digits indicate nozzle area ratio. 
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thrustors are identified in table I by the first digit of the unit designation. The number 3 
refers to a unit with a 1/2-inch-diameter heater coil with heater supports made of boron 
nitride. The number 4 corresponds to a 1-inch-diameter heater coil with heater supports 
made of tungsten. The letters W and BN refer to tungsten or boron nitride nozzle 
bodies, respectively. Only the 3W25 unit has a tungsten nozzle body. All other nozzle 
bodies were made of boron nitride to simplify fabrication. The boron nitride nozzle 
bodies performed satisfactorily for the comparatively short runs and moderate tempera- 
tures required to obtain the necessary performance data. The remaining digits indicate 
the nozzle area ratio. 


Range of Test Parameters 

Data were obtained on each of the seven different nozzles at propellant flow rates of 
-4 -4 -4 

0. 5x10 , 1. 0x10 , and 2. 5x10 pound per second. For each of these propellant flow 

rates, specific power (input electrical power /propellant flow rate) values of 0, 0. 5x10^, 

1 . 0x10^, and 1 . 5x10^ kilowatts per pound per second were used except for the lowest 
propellant flow rate of 0. 5x10”^ pound per second. At this flow rate, a single specific 
power point of approximately 1. 2x10^ kilowatts per pound per second was used corre- 
sponding to the minimum current output of the power supply for the heater coils used. No 
higher specific power points were attempted because of the marginal thrustor cooling at 
this low flow rate. 

4 

The maximum specific power value was limited to 1. 5x10 kilowatts per pound per 
second to minimize the component problems that are encountered at the higher propellant 
temperatures. Each specific power point was held for a period of time in excess of 
15 minutes, during which time the measured stagnation pressure was monitored to assure 
that equilibrium conditions had been achieved. 


Research Facility and Instrumentation 

The experimental program was conducted in one of the large vacuum tanks at the 
Lewis Research Center. This facility (described in detail in ref. 12) has a diameter of 
15 feet and a length of 63 feet. Exhaust pumping speed capabilities of 12 000 cubic feet 
per minute provided ambient pressures below 0. 3 millimeter of mercury for all pro- 
pellant flow rates. The low ambient pressures were more than sufficient to ensure com- 
plete expansion for all the nozzle area ratios studied. Electrical power was supplied by 
a commercial solid-state conversion d-c power supply. Gaseous hydrogen was supplied 
from a bank of four commercial hydrogen bottles. 
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Figure 6. - Schematic drawing of thrust stand and instrumentation. 


A schematic drawing of the thrustor 
mounting and system instrumentation is 
shown in figure 6. The propellant flow 
rate was metered by a calibrated choked 
orifice with strain-gage pressure trans- 
ducers to measure orifice upstream and 
downstream pressures. Thrustor cham- 
ber pressure was measured with a 
Bourdon tube-type absolute-pressure 
gage and a strain-gage pressure trans- 
ducer. Voltage and current measure- 
ments were made with conventional in- 
strumentation. Thrust measurements 
were made with a four-point suspension- 
platform-type thrust stand (ref. 13). A 
differential transformer immersed in an 
oil bath, to maintain a constant- 
temperature environment, was used as 
the thrust-sensing element. 

The overall measurement accuracy 
of the data is as follows: chamber pres- 


sure and propellant flow rate measure- 
ments are accurate to within ±2 percent. Thrust measurements as determined by weight 
calibration before each test run under vacuum conditions are accurate to within ±1 milli- 


pound. The overall accuracy of the measured thrust coefficient is estimated to be ap- 
proximately ±10 percent at the low-thrust levels and ±3 percent at the high-thrust levels. 


DATA REDUCTION 

The measured data, along with other computed values, are presented in tables I 
and II for each of the seven different nozzle configurations. Table I (p. 8) presents the 
cold-flow data (obtained without heat addition) and table n presents the hot-flow data (ob- 
tained with heat addition to the propellant). 

The measured thrust coefficient and the corresponding loss in thrust coefficient are 
calculated as 


C 


F “ 


m 


P m A th 


( 9 ) 
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TABLE n. - HOT- FLOW DATA 


Unit 

designa- 

tion 

(a) 

Nozzle 

area 

ratio, 

V A th 

Nozzle 

throat 

flow 

area, 

A th’ 
sq in. 

Isentropic 

thrust 

coefficient, 

° F i 

Propellant 
flow rate, 
w, 

lb/sec 

Measured 

total 

pressure, 

P m- 

psia 

Measured 

thrust, 

F m’ 

lb 

Measured 

thrust 

coefficient, 

C F 

m 

Difference 
between 
calculated 
isentropic 
and mea- 
sured 
thrust 
coefficient, 

Cf * 

Throat 

Reynolds 

number, 

Re th 

3W25 

25 

3. 63x10" 3 

1. 652 

2. 5xl0" 4 

21. 5 

114. 5x10" 3 

1.468 

0. 184 

4770 





2. 5 

23.4 

124. 2 

1.463 

. 189 

4270 





1. 0 

9.4 

43.4 

1. 280 

. 372 

1785 





1.0 

10.4 

47.8 

1. 267 

. 385 

1565 





1.0 

9.0 

42.7 

1. 308 

. 344 

1895 





1.0 

10. 35 

48.9 

1.300 

. 352 

1575 





1.0 

10. 5 

48. 7 

1. 278 

. 374 

1550 





1.0 

11. 3 

51.6 

1.259 

. 393 

1400 





. 5 

4.9 

18.7 

1.051 

. 601 

867 





. 5 

5.0 

19.6 

1.080 

. 572 

844 





. 5 

5. 1 

19. 8 

1.070 

. 582 

823 





. 5 

4.9 

18.05 

1.014 

. 638 

867 

3 BN 50 

50 

3. 84x10" 3 

1.677 

2. 5X10" 4 

20. 6 

113. 0xl0 -3 

1.430 

0. 247 

4740 





2. 5 

25. 3 

137.0 

1.410 

. 267 

3610 





2. 5 

28. 5 

153.0 

1.398 

. 279 

3070 





1.0 

8. 75 

42. 3 

1. 260 

.417 

1820 





1.0 

10. 3 

48.9 

1. 235 

.442 

1465 





1.0 

11.05 

51.8 

1. 220 

.457 

1330 





1.0 

11. 75 

53. 8 

1. 193 

.484 

1225 





. 5 

4.4 

18.05 

1.068 

. 609 

903 

3BN75 

75 

3. 63x10" 3 

1.691 

2. 5xl0" 4 

22.0 : 

118. 5xlO" 3 

1.485 

0. 206 

4750 





2. 5 

25. 7 

136. 5 

1.462 

. 229 

3890 





2. 5 

28. 9 

153.0 

1.458 

. 233 

3300 





1. 0 

9.4 

46. 1 

1. 350 

. 341 

1780 





1. 0 

11.0 

50.3 

1. 261 

.430 

1445 





1. 0 

11. 8 

52. 7 

1. 230 

.461 

1330 





. 5 

4.9 

18. 5 

1.040 

. 651 

890 

3BN100 

100 

4. 07x10" 3 

1.696 

2. 5xl0" 4 

21.4 

123. 6X10" 3 

1.420 

0. 276 

3940 





2. 5 

24. 8 

142.6 

1.413 

. 283 

3220 





2. 5 

25.9 

145. 5 

1. 381 

. 315 

3050 





2. 5 

20. 5 

121.3 

1.455 

. 241 

4170 





2. 5 

24.9 

141.4 

1. 396 

. 300 

3200 





1.0 

8.4 

41.9 

1. 225 

.471 

1755 





1 . 0 

9. 82 

50.8 

1.272 

.424 

1430 





1 . 0 

11. 2 

54. 2 

1. 190 

. 506 

1195 





1.0 

9.9 

47.9 

1. 190 

. 506 

1405 





1.0 

10. 15 

53.3 

1. 173 

. 523 

1205 





. 5 

4. 35 

17.43 

.985 

. 711 

880 





. 5 

4.40 

18.05 

1.010 

. 686 

867 


a First digit: 3 has 1/2-in. -diam heater coil and boron nitride heater supports, while 4 has 1-in. -diam heater 
coil and tungsten heater supports; letters W and BN indicate tungsten or boron nitride nozzle body; last 
digits indicate nozzle area ratio. 



TABLE H. - Concluded. HOT- FLOW DATA 


Unit 

designa- 

tion 

Nozzle 

area 

ratio, 

V A th 

Nozzle 

throat 

flow 

area, 

A th’ 
sq in. 

Isentropic 

thrust 

coefficient, 

Propellant 
flow rate, 
w, 

lb/sec 

(a) 





4BN100 

100 

4. 65X10' 3 

1. 676 

2. 5xl0" 4 
2. 5 
2. 5 
1.0 
1.0 
1.0 
1 . 0 
. 5 
. 5 

4BN40 

40 

11. 90X10' 3 

1. 670 

2. 5X10" 4 
! 2. 5 
2. 5 
1.0 
1.0 
1.0 
. 5 
. 5 

3BN150 

150 

1. 32x10* 3 

1 

1. 705 

1 

2. 5X10" 4 
2. 5 
2. 5 
2. 5 
1.0 
1.0 
1.0 
. 5 


Measured 

Measured 

Measured 

Difference 

Throat 

total 

thrust, 

thrust 

between 

Reynolds 

pressure, 

F m> 

coefficient, 

calculated 

number, 

p m> 

psia 

lb 

C F 

m 

isentropic 
and mea- 
sured 
thrust 
coefficient, 

C F 

r S. 

Re th 

15.9 

113. 4xl0" 3 

1. 534 

0. 162 

4310 

21. 1 

147.0 

1.497 

. 199 

2970 

23. 3 

160.4 

1. 480 

. 216 

2610 

8. 5 

48. 7 

1. 232 

.464 

1235 

7.9 

48. 1 

1. 310 

. 386 

1370 

7.2 

45. 2 

1. 350 

. 346 

1545 

9. 1 

54.4 

1. 285 

. 411 

1128 

3.4 

18. 3 

1. 157 

. 539 

895 

3. 8 

19. 6 

1.110 

I 

. 586 

773 

6.6 

110. 7xl0” 3 

1.408 

0. 262 

2650 

8. 35 

135. 2 

1. 360 

. 310 

1945 

9.6 

154.0 

1. 348 

. 322 

1615 

2. 70 

39. 1 

1. 217 

. 453 

1095 

3.35 

46. 6 

1. 170 

. 553 

838 

3.90 

52. 3 

1. 127 

. 543 

680 

1.70 

18. 54 

. 917 

. 753 

443 

1.78 

18. 95 

. 895 

. 775 

425 

59.4 

118. 5x10” 3 

1. 512 

0. 193 

7740 

59.05 

116. 8 

1.496 

. 209 

7800 

73. 8 

145.4 

1. 491 

. 214 

5780 

78.4 

158.4 

1. 532 

. 173 

5360 

24.4 

43.0 

1. 335 

. 370 

3050 

28. 6 

48. 5 

1. 285 

. 420 

2460 

31. 85 

53. 0 

1. 260 

. 445 

2140 

13.0 

17. 86 

1.041 

. 664 

1445 


a First digit: 3 has 1/2-in. -diam heater coil and boron nitride heater supports, while 4 has 1-in. -diam heater 
coil and tungsten heater supports; letters W and BN indicate tungsten or boron nitride nozzle body; last 
digits indicate nozzle area ratio. 
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( 10 ) 


° F t = Cf 1 ' CF m 


The throat Reynolds number is defined as 


Re th " 


w d th 
A th^th 


( 11 ) 


Since the propellant temperature was not measured (except at the cold-flow conditions), 
the viscosity for the hot-flow data could not be explicitly evaluated. Therefore the follow- 
ing approximate method was used to determine the throat Reynolds numbers. 

For the cold-flow data, the viscosity was evaluated at the throat temperature cor- 
responding to an isentropic process from the known stagnation temperature of 530° R. 

For the hot data, the throat Reynolds number ratio (Re^ co id/R e th hoP R>r same 
nozzle throat diameter and the same propellant flow rate is given by 


Re th, hot _ ^th, cold 
Re th, cold ^th, hot 


(12) 


For the temperature range covered, an adequate expression for the viscosity of hy- 
drogen (fitted to the data of ref. 3) is /i oc Equation (12) can then be expressed 

as 


Re th, hot 
Re th, cold 


T \2/3 
i th, cold \ 

, T th, hot / 


(13) 


For a one-dimensional isentropic process, the temperature at the nozzle throat is 
proportional to (P^A^/w) • Ror conditions in which the nozzle throat area and propel- 
lant flow rate are constant, the temperature ratio in equation (13) is then given by 


T th, cold 
^th, hot 



(14) 


If it is assumed that the pressure at the nozzle throat and the measured chamber pressure 
are similarly related for conditions of both hot and cold flow as 
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the ratio of the cold to the hot throat Reynolds number can be approximated by 

Re th, hot ^ Am, cold 
Re th, cold \ P m, hot 

Although this approximate procedure for estimating the propellant temperature may in- 
troduce inaccuracies, the resulting error is considered to be within the experimental 
accuracies of the other parameters. 



EXPERIMENTAL RESULTS 

The experimental data are suitable for the determination of thrustor performance 
but they are not suitable for a precise interpretation of the nature or the mechanism of 
the nozzle losses. This statement is especially applicable to the hot data, which includes 
both viscous and heat-transfer effects, that cannot be separated from each other. The 
data are useful, however, to identify correlating parameters for analytically predicting 
overall performance. 

Figure 7 presents the thrust coefficient for the cold and the hot data for the 3W25 
nozzle. The trends exhibited by these data are typical of all the other nozzles. The de- 
creasing thrust coefficient with decreasing throat Reynolds number is consistent with the 
predictions of the laminar-flat-plate analysis, but the magnitude of the measured thrust 
coefficients is much less. This was not unexpected and indicates that the conical nozzle 
losses are much higher than the calculated laminar-flat-plate losses. The magnitude of 
the measured thrust coefficients at the low throat Reynolds numbers is also significant. 
The measured thrust coefficient decreases from a value of approximately 1. 45 at 
Re th w to a va * ue approximately 1. 05 at Re th « 850. This latter value is less 

than the calculated isentropic thrust coefficient (1. 27) for a choked orifice. 

In figure 7, the hot-flow thrust coefficients are higher than the cold-flow values at a 
given throat Reynolds number. If equation (7) is used as a guide to estimate the differ- 
ence between the hot and cold data, the experimental trend is contrary to that predicted 
by the analysis. The only parameter in equation (7) that changes as a result of the 
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Figure 7. - Typical comparison of hot and cold thrust coefficients. Unit 
designation, 3W25. 


differences between the hot and cold operating conditions is the factor f , which depends 

on the wall temperature. For cold flow, the expansion process is essentially adiabatic, 

and therefore T w /T q ~ 1. 0. For hot flow where cooling of the nozzle walls occurs by 

radiation, T /T < 1. 0. As T /T decreases, the loss in thrust coefficient, as given 
w o w o 

by equation (7) (including viscosity effects), increases and correspondingly the thrust 
coefficient for the hot-flow data at a given throat Reynolds number should be less than 
that for the cold data. The experimental data therefore indicate that the throat Reynolds 
number based on the local stream temperature is not a good correlating parameter both 
for the hot- and cold-flow data. 

Even though the free- stream throat Reynolds number would not correlate both the hot 
and cold data, it seemed probable that a suitable correlation could be developed for the 
loss in thrust coefficient if the throat Reynolds number were based on an appropriate ref- 
erence temperature. The nozzle-wall temperature was found to be an acceptable choice. 
The throat Reynolds number based on the wall temperature is given by 


Re 


th, w 



(16) 


In terms of the previous throat Reynolds number, equation (16) becomes 


Re 


th, w 




(17) 


For the conventional boundary -layer assumption of zero radial pressure gradient and 


a viscosity-temperature relation for hydrogen of p oc T 


2/3 


Re th, w is given by 
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(18) 


Re 


th, w 



For the cold-flow data, the nozzle-wall temperature can be approximated by the adia- 
batic wall temperature. Thus, for a laminar recovery factor of 0. 84, 


I— J = 0. 857 (19) 

' w 'th, cold 

For the hot data, an estimate of (T/T ) is not as straightforward. Optical pyrom- 

th 

eter measurements indicated that the nozzle-wall temperature for similar operating con- 
ditions was of the order of 0. 50 < T w /T e < 0. 66, where T E is the ideal temperature 
the propellant would attain in the heater chamber for a heating efficiency of unity (i. e. , 
if all the input electrical power were converted to thermal energy). If T-, is introduced, 

ili 



th, hot 



( 20 ) 


If the value of T^/T^, is taken as the one-dimensional isentropic value of 0. 833 and 
T w, tt/ T E is approximated as th /T E » T w /T e » 0. 6, an estimate for (T/T w ) 



0. 833 
0. 60 


1. 388 


( 21 ) 


Although a constant value of (T/T ) is unlikely for the complete range of op- 

w th, hot 

erating conditions covered herein, a more accurate value could not be obtained in this 
investigation. 

Figure 8 is a comparison of the hot and cold losses in thrust coefficient in terms of 

Re th w * or each the seven different nozzles. The data indicate, in general, that, at 

the high propellant flow rates where the heat-transfer losses and thermal -boundary layer 

are small, the hot and cold loss in thrust coefficient are of the same magnitude. As the 

propellant flow rate decreases, the losses for the hot data tend to increase over that of 

the cold-flow data at the same value of Re.. 

th, w 
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(b) Nozzle area ratio, 50, unit designa- 
tion, 3BN50. 



(c) Nozzle area ratio, 75; unit designa- 
tion, 3BN75. 
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(d) Nozzle area ratio, 100, unit designa- (e) Nozzle area ratio, 100; unit designa- 
tion, 3BN100. tion, 4BN100. 
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Wall temperature throat Reynolds number, Re^ 

(f) Nozzle area ratio, 40, unit designation, 4BN40. (g) Nozzle area ratio, 150, unit designation, 3BN150. 


Figure 8. - Comparison of hot and cold data in terms of wall temperature throat Reynolds number. 





Figure 9. - Correlation of loss In thrust coefficient for cold-flow data. 


The curves shown in figure 8 have the functional form Cp oc ( Re th 1//2 a 1 * 1 are 

in, w 

a consistent and reasonable approximation of the data, except for the 4BN100 nozzle. The 
data obtained with the 4BN100 nozzle were much different from the data obtained with the 
3BN100 nozzle. In addition, the slope of the line to fit the 4BN100 nozzle data was not 
consistent with the data from all the other nozzles. No explanation can be given for this 
difference. 

All the data (except that for the 4BN100 nozzle) presented in figure 8 exhibited a con- 
sistent trend of increasing losses with increasing nozzle area ratio. The area-ratio 
correlation can best be shown in terms of the cold-flow data as presented in figure 9. 
Based on a straight-line correlation of the data, the final relation for the loss in thrust 
coefficient is given by 


C 



17. 6 exp 0. 0032 



( Re th, w^ 


1/2 


The effectiveness of this correlation to estimate the thrust coefficient for both cold- 
and hot-flow conditions for the area ratios investigated is shown in figure 10. The data 
lie in a ±5 -percent band, except for the data from the 4BN100 nozzle and the data at the 
lowest propellant flow rate. It is believed that the reason the data at the lowest propel- 
lant flow rate deviate from the correlation is that the heat -transfer losses become more 
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important at low throat Reynolds num- 
bers. However, these data at the low 
propellant flow rates should not be 
considered as general but rather as 
characteristic of the radiation-cooled 
thrustor used in this study. Heat- 
transfer losses are governed by de- 
sign considerations. Other nozzle- 
cooling methods, sizes, and mate- 
rials could yield different results. 

The approximations made to ob- 
tain the final correlation of the data 
may not be suitable for other ranges 
of propellant flow rates and specific 
power levels than those investigated 

.8 .9 1.0 i.i 1.2 1.3 1.4 1.5 1.6 herein. For example, it has been 

Calculated thrust coefficient from correlation . 

assumed that T /T-, = 0. 6 for the 

Figure 10. - Comparison of experimental thrust coefficients with values w & 

calculated from correlation. hot data and that T ,/T™ = 1. 0 for 

the cold-flow test data. This temper- 
ature ratio, of course, must vary continuously with the propellant temperature. Further- 
more, the nozzle-wall temperature has an axial variation and is not constant as was as- 
sumed. Additional research is required to refine the correlation and to understand better 
the various mechanisms contributing to the losses in low-thrust-level nozzles. 

CONCLUSIONS 

The following conclusions regarding the performance losses of nozzles at low-thrust 
levels have been drawn from the experimental data: 

1. Nozzle losses can seriously degrade the thrust coefficient of nozzles with low 
throat Reynolds numbers to values much below those calculated for an isentropic expansion 
process. In fact, thrust coefficients were measured that were less than the calculated 
thrust coefficient for isentropic flow through a choked orifice. 

2. The measured nozzle losses are significantly greater than the viscous losses cal- 
culated from the skin friction relation for laminar flow over a flat plate. The trends, 
however, in terms of increasing losses with decreasing throat Reynolds numbers are 
similar. 
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3. A qualitative correlation of the loss in thrust coefficient for both the hot and cold 
data was obtained in terms of the nozzle area ratio and throat Reynolds number in which 
the temperature-dependent properties (density and viscosity) are evaluated at the nozzle 
wall temperature. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, July 26, 1965. 
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APPENDIX A 


SYMBOLS 


A 

nozzle flow area 

T 

A / A th 

geometric nozzle area ratio 

t e 

( A / A th ) 
tn eff 

effective nozzle area ratio 


C 

factor of proportionality in 



eq. (B7), njn = C(T W /T) 

u 

C F 

thrust coefficient, F/P q A^ 

w 

C F, 

difference between calcu- 

X 

a 

lated isentropic thrust 
coefficient and measured 
thrust coefficient 

a 

y 



6* 

% 

loss in thrust coefficient 
due to boundary-layer vis- 

A 


cous losses, F v /P Q Ay 1 


d 

nozzle diameter 

M 

F 

thrust 

P 

F i 

loss in thrust due to 

r 


nozzle losses 

Subscripts; 

F 

V 

loss in thrust due to vis- 

e 


cous drag on nozzle sur- 
face 

ht 

g 

gravity acceleration 

i 

M 

Mach number 

m 

P 

total pressure 

o 

P 

pressure 


Re 

Reynolds number 

th 

r 

nozzle radius 

V 

S 

nozzle surface 

w 

s 

distance along nozzle 

X 


surface 



temperature 

ideal temperature attained 
by propellant for heating 
efficiency of unity 

fluid velocity 
propellant flow rate 
axial distance along nozzle 
nozzle divergence angle 
specific heat ratio 
displacement thickness 

nozzle angle divergence coef- 
ficient, A = (1 + cos &)/2 

fluid viscosity 

fluid density 

shear stress 

nozzle-exit conditions 
heat-transfer effects 
isentropic conditions 

experimentally measured 
values 

stagnation conditions up- 
stream of nozzle throat 

throat conditions 

viscous effects 

wall conditions 

axial component 
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APPENDIX B 


CONICAL NOZZLE VISCOUS LOSSES 


It will be assumed in the following development that (1) the important viscous losses 
take place in the supersonic portion of the nozzle, (2) the boundary layer begins at the 
nozzle throat, (3) the nozzle is an isothermal surface, and (4) the local flow conditions 
can be obtained in terms of the stagnation conditions from one-dimensional isentropic- 
flow relations. The divergent portion of a conical nozzle with pertinent coordinates is 
shown in figure 1 (p. 3). 

The loss in thrust due to the axial component of the viscous drag is given by the sur- 
face integral of the local shear stress over the nozzle surface: 




2nrr cos a ds 


(Bl) 


For purposes of calculation, it is desirable to express equation (Bl) in terms of the 
geometric nozzle area ratio A/A^. The linear coordinates are related to the nozzle 
area ratio as follows: 


x 



1 


d th 2 tan a 


(B2) 



= ds cos a 


(B3) 


Equation (Bl) may now be written in the form 

f A e /A th 

The loss in thrust can be expressed in terms of a loss in thrust coefficient C„ by 

l 1 ,. 


tA. 


tan a 


^ d(— 


Hh/ 


(B4) 
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I 


c = i. 

F v P.A 


oth 


P Q tan a 



(B5) 


The loss in thrust coefficient can be computed by substituting into equation (B5) the ap- 
propriate expression for the local shear stress r and then performing the indicated in- 
tegration. Unfortunately no useful relations for the laminar shear stress on a nozzle sur- 
face are known. However, the skin friction relation for laminar flow over an iso- 
thermal flat plate provides a simple initial approximation for comparison with exper- 
iment. This expression, as given in reference 9, is 



0. 664 C°- 5 



(B6) 


The constant C relates the viscosity and temperature at the wall and free- stream con- 
ditions: 


V 



(B7) 


The local length Reynolds number in equation (B6) can be expressed in terms of the noz- 
zle throat Reynolds number Re^ as follows: 


— = Re a,- 6 -— — — (B8) 

p p th u th d th 11 


By suitable substitution of the one-dimensional isentropic-flow relations and equation (B2) 
into equation (B8), the local length Reynolds number can be written as 
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y 

Y - 1 


(BIO) 


- pu 2 = 1 PM 
2 2 ° 


2 



Substituting equations (B7), (B9), and (BIO) into equation (B6), and the resulting ex- 
pression for the shear stress into equation (B5) yields 


C 

F 

v 


0. 664 f 


(Re th tan a) 


fcl. 


V A th 


f l f 2 d 


th/ 


(Bll) 


where 



The local stream properties M and T/T in the parameter f 0 are evaluated from one- 

o z 

dimensional isentropic relations for the local geometric nozzle area ratio (uncoupled case) 
or at an effective nozzle area ratio (coupled solution), which includes the influence of the 
boundary- layer growth on the free-stream conditions. 

For the coupled solution, the effective nozzle area ratio is defined in terms of the 
boundary-layer displacement thickness and geometric nozzle area ratio by 



(B12) 
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where 



(B13) 


For laminar flow over a flat plate with constant wall temperature, 6 */ x may be de- 
rived from reference 9 as 

= c 1 / 2 1//2 [l. 73 + 1. ll(r- 1)M 2 ] (B14) 

where C and the local length Reynolds number are given by equations (B7) and (B9). 
Hence, 5*/r may be expressed as 


5 * _ goglfg 
r (Re th tan Q!) 1//2 


(B15) 


where 
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Equation (B15) must be solved by iteration to account for the change in free-stream 

conditions in gg due to the boundary- layer growth. As a first approximation, equation 

(B15) may be solved assuming (A/A,. ) = A/A,. . This gives a value of 6*/ r that per- 

tn eff tn 

mits equation (B12) to be solved for (A/A^) , which becomes the basis for the next cal- 

eff 

dilation. If the boundary layer is thin, this process converges rapidly but for thick 
boundary layers the process is facilitated by using a graphical solution. 
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